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We present time-series radial-velocity and photometric observations of a solar- 
type double-lined eclipsing binary star (V12) in the old open cluster NGC188. 
We use these data to determine the spectroscopic orbit and the photometric 
elements for V12. From our analysis we determine accurate masses {Mp = 
1.103 ± 0.007 Mq, Ms = 1.081 ± 0.007 M©) and radii (Rp = 1.424 ± 0.019 Rq, 
Rs = 1.373±0.019 Rq) for the primary (p) and secondary (s) binary components. 
We adopt a reddening of Eb-v = 0.087 for NGC 188, and derive component ef- 
fective temperatures of 5900 ± 100 K and 5875 ± 100 K, respectively, for the pri- 
mary and secondary stars. From their absolute dimensions, the two components 
of V 12 yield identical distance moduli of Vq — My = 11™24±0^09, corresponding 
to 1770 ± 75 pc. Both stars are near the end of their main-sequence evolutionary 
phase, and are located at the cluster turnoff in the color-magnitude diagram. We 
determine an age of 6.2 ± 0.2 Gyr for V 12 and NGC 188, from a comparison with 
theoretical isochrones in the mass-radius diagram. This age is independent of 
distance, reddening, and color-temperature transformations. We use isochrones 
from Victoria-Regina and Yonsei-Yale with [Fe/H] = -0.1 and [Fe/H] = 0.0. From 
the solar metallicity isochrones, an age of 6.4 Gyr provides the best fit to the bi- 
nary components for both sets of models. For the isochrones with [Fe/H] = 
-0.1, ages of 6.0 Gyr and 5.9 Gyr provide the best fits for the Victoria-Regina and 
Yonsei-Yale models, respectively. We use the distance and age estimates for V 12, 
together with best estimates for the metallicity and reddening of NGC 188, to 
investigate the locations of the corresponding VRSS and Y"^ isochrones relative 
to cluster members in the color-magnitude diagram. Plausible changes in model 
metallicity and distance to better match the isochrones to the cluster sequences, 
result in a range of ages for NGC 188 that is more than 3 times that resulting 
from our analysis of V 12. 

Subject headings: Open clusters: general — Open clusters: individual NGC 188 

— Stars: binaries: spectroscopic — Stars: binaries: eclipsing — Stars: evolution 

— Techniques: spectroscopy — Techniques: photometry 



Introduction 



NGC 188, once thought to be the oldest open cluster in our Galaxy (e.g.. lSandagelll962l ). 
represents the old stellar population of the Galactic disk. Due to its richness, location in 
the Galaxy (02000 = 0^ 47 "^, S^mn = +85° 15^ ; / = 122.8°, 6 = 22.4°), large population of 



single and binary members (jPlatais et al. II2003I : iGeller et al. II2008I ). and its age, NGC 188 is a 



- 3 - 



benchmark cluster for studies of stellar evolution and cluster dynamics, and for the formation 
and chemical and dynamical evolution of the Galactic disk. Determining an accurate and 
precise age for NGC 188 is of general importance. 

The age determined for NGC 188 has changed dramatically si nce the early determination 



by Sandage ( 1962 ) of 14-16 Gyr. Later studies derived 9- 12 Gyr (Demargue fc Larson! 



Iben 



19671'). and more recent values range from ~6 Gyr 



Caputo et al. Ill990l: iDemargue et al. 



1992 



Mevnet et al. 



Twarog fc Anthonv-Twarog 



1993 



Dinescu et al. 



IVandenBerg fc StetsonI I2OO4I ). Past 



1964 



1989 



9951) to ~7- 



8 Gyr flHobbs et al. 



1990 



Sarajedini et al. 



1999 



age 

estimates were based on fitting model isochrones to the cluster's main sequence, turnoff, and 
giant branch in the color- magnitude diagram (CMD). Ages derived using this "isochrone 
method" are sensitive to uncertainties in cluster extinction and reddening, as well as in the 
conversions between model (bolometric luminosity, temperature) and observed (apparent 
luminosity, color) quantities. 

For NGC 188, the values for reddenin g and metallicity are w ell constrained. A redden- 
ing of E{B — V^) = 0?^087 der i ved fr om the ISchlegel et al. I (119981 ) dust maps agrees with the 
estimates of ISarajedini et al.l (119991 ) of 0™09 ± 0™02 from their two-color diagram. We use 
E{B -V) = 0^087 throughout t his paper. The m e talhcity of NGC 188 has been studied 
during the last two decades (e.g. iHobbs et al.l Il990l : iFriel &: Janes I Il993l : iFriel et al. I l2002l : 
' 2003h . Results from photometric studies agree with 



Randich et al. II2003I : IWorthev fc Jowett 



results from medium- to high-resolution spectroscopy, both of which find the cluster's metal- 
licity to be in the range from [Fe/H] = —0.12 to solar. Estimates of the cluster distance 
from isochrone-fitting to the CMD, ha ve over the past decade agreed to within a few per 



cent, ranging from 1660 pc to 1710 pc (Sarajedini et al. Ill999l : IVandenBerg fc StetsonI 12004 
Bonatto et al. 2005: Fornal et al. 20071) with individual uncertainties of the order of 5%. 



NGC 188 is a key cluster of the VV IYnEI Open Cluster Study (WOCS; iMathieu I [200o[ ) . 
As part of WOCS, iGeller et al.l (120081 ) present a time-series radial-velocity survey of more 
than 1046 stars in the field of NGC 188 ov er a time baseline of 11 years. In combination with a 
deep proper- motion study of the cluster (jPlatais et al.ll2003l ). the WOCS spectroscopic data 
have securely identified large populations of single and binary cluster members. V 12 was 
established as a double-lined spectroscopic binary and confirmed as a clu ster member early in 
the WOCS survey. The eclipsing nature of the system was discovered bvlZhang et al.l (120021 ) 
in a time-series photometric survey for variable stars. IZhang et al.l (120021 ) recorded only one 



^The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, 
Yale University, and the National Optical Astronomy Observatories located at Kitt Peak National Observa- 
toty, Arizona, USA. 
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co mplete eclipse and c lassified V 12 as an Algol- type eclipsing binary. The photometry used 
bv Izhang et aP J2002h locate d V12 at the foot o f the giant branch in the NGC188 CMD. 
A later photometric survey by IZhang et al. I (120041 ) placed V12 above the cluster turn-off in 
the CMD, and estimated the period of V 12 to be 2.8 days. Just recently , V 12 was identified 
as a detached eclipsing binary in NGC 188 by iMochejska et al.l (120081 ). They determine a 
period of 2.9 days. 

As a spectroscopic double-lined and eclipsing binary, precise masses and radii can be 
determined for the primary and secondary components of V12. Because both components 
are near the turnoff mass for NGC 188, a comparison of their masses and radii to theoretical 
isochrones can put tight constraints on the age of the binary, and thus its parent cluster. 
Importantly, the determinations of the masses and radii are independent of the reddening, 
distance, and metallicity of the system. In addition, double - lined eclipsing binaries are 
excellent primary distance indicators (e.g. iGuinan et al. Ill998l : IClausen II2004I ). Therefore, 
the analysis of eclipsing binary stars, in particular when coupled with reliable information 
about the stellar chemical composition, is important to improving our understanding of the 
physical properties and evolution of stars. However, stellar evolution models have progressed 
to a level where stellar parame ters must be obs erved with a precision of ~1% or better to 
provide useful constraints (e.g. I Andersen Ill99ll ). For examples of det ailed analysis of fiel d 

J2006h . 



binaries and the comparisons with stel lar ev olution models see e. g 
Lopez- Morales. M. fc Shaw. J. S.I (12003), and lciausen et al.l koQ^ . 



Torres et al. 



While bright eclipsing binaries in re 


atively nearby clusters have been observed for 


decades (Kraft & Landolt 


1959; 


Semeniuk 


1962 


Sahade & Davila Ill963: Popova & Kraicheva 


1984: Clausen & Gimenez 


1987 


Gimenez & C 


ausen 


1996] 


, and several studies have been 


published recently (see e.g 


Southworth et al. 


( 


2005 


) and 


Southworth & Clausen (2006) 



and references therein), it is only re cently that fainter sy s tems in more distant clusters have 
been discovered and published (e.g. iKaluzny et al.ll2006l : iGrundahl et al.ll2008l ). 



In Section [2] we begin the analysis of the spectroscopic data on V12. We determine 
the orbital elements of V 12 and give the probabilities for it being a kinematic member of 
NGC 188. Section [3] describes the photometric data, the light curves, and the photometric 
elements of the V 12 components. That leads to our determination of the physical parameters 
for V12, including effective temperature, distance, and age, in Section HI In Section [5] we 
examine how well the isochrones that best fit the V 12 components represent the rest of the 
cluster members in the NGC 188 CMD, and we compare the cluster age determined via V 12 
to those determined by the isochrone method. In Section we summarize our results and 
give our conclusions. 
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2. Radial velocities and spectroscopic elements of V 12 
2.1. Spectroscopic Observations and Datareduction 



A detailed description of the spectrosc opic observat i ons a nd radial-velocity measure 



ments of stars in NGC188 can be found in iGeller et al. I (120081 ) . The spectroscopic obser- 



vations of V 12 were acquired from Sep 1996 to July 2001 with the WlYN telescope and 
the Hydra multi-object spectrograph in a single 250A wide echelle order centered on 5130 
A and with a resolution of ~ 20 000. For each observation three 40 minute exposures were 
combined to increase the signal-to-noise ratio of the spectrum and to filter out cosmic ray 
events. The CCD bias level was determined using the overscan-section of the chip and 
subtracted. A single domefiat exposure, adjacent in time to the science exposures, was 
used for the purposes of tracing and extracting the stellar spectrum, fiatfield correction, 
and fiber throughput correction. The three consecutive science exposures were bracketed by 
Thorium- Argon (ThAr) comparison lamp exposures used to define the pixel-to-wavelength 
mapping (dispersion solution). The DOHYDRaI^ IRAF task was used for the reduction of 
the WIYN/Hydra multi-spectrum data. The signal-to- noise ratio (S/N) for the 16 spectra 
used for the orbital solution of V12 ranges from ~15-25 per resolution element (2 pixels) of 
~ 14 kms~^. For single lined stars wi th slow to rnodera te rotation velocities, radial veloci- 
ties with a precision of < 0.4 kms~^ (IGeller et al. II2008I ) were derived from the spectra via 
cross-correlation with a high S/N sky (solar) spectrum obtained at evening twillight. 



2.2. Component Radial- Velocities and Orbital Elements 



Because we seek maximum precision and accuracy for the V 12 velocities, we have used 
two methods to independently re-reduce the spectra and derive velocities for both compo- 
nents. First, we have determined radial veloci ties using the two-dimensional cross-correlation 
algorithm TODCOR JZucker fc Mazeh Ill994h . For the TODCOR analysis we have used iden- 
tical synthetic templates for the two binary components, since they are nearly identical; see 
Section m The templates were calculated from the Synplot/Syn^peB tool by Hubeny & Lanz 
for Teff = 6000 K, log(^) = 4.0, [Fe/H] = 0.0, ^;s ini = 20 km s-\ and a resolution of R = 
20 000, adopting an ATLAS9 model atmosphere ( Kurucz1ll993l ). Preliminary spectroscopic 
elements were calculated from the measured radial veloc ities using the method of Lehman- 
Filhes implemented in the SBOP program (lEtzel 1120041 ). The orbital period was adopted 



^http:/ /iraf.noao.cdu/tutorials/dohydra/dohydra.html 
■^http: / /nova. astro. umd.edu/Synspec43/synspec. html 
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from Equation [2] (Section 13.31) . The orbit was assumed to be circular based on an initial ec- 
centricity of 0.005 ±0.003. Because spectral lines of the stars move in and out of the observed 
spectral range with orbital phase, it is nece ssary to apply corre ctio ns to the radial velociti es 
derived from the TODCOR algorithm (e.g. [Torres et al. I (119971 ) and [Torres &: Ribad (120021 )). 
To determine such corrections, we have constructed synthetic binary spectra from the tem- 
plate, adopting for each orbital phase the radial-velocity shifts defined by the preliminary 
elements. These synthetic spectra were then analyzed by TODCOR in the same manner as 
the observed V 12 spectra. The corrections (i.e. measured — calculated velocities for the 
synthetic binary spectra), which are shown in Figure [T|, are then added to the measured 
V12 velocities and revised spectroscopic elements are determined from SBOP. In the case 
of V12, we find that the radial-velocity corrections are quite significant, up to 3 kms~^ for 
both components. Including the corrections improves the fit of the theoretical orbit to the 
radial velocities significantly, and the effect on the absolute masses is not negligible (~2% for 
both components). The final velocities including corrections are listed in Tabled and the 
corresponding spectroscopic elements are given in Table [2l The final spectroscopic orbital 
solution for V 12 is shown in Figure [2l Within the errors, separate analysis of the radial 
velocities for the primary and secondary components, respectively, yield identical elements. 
Also, analysis performed with several other templates (f sini = 10.0 kms~^, correspond- 
ing t o synchronous rotation (see Section 1^^ : 200 K lower tempe ratures; different spectrum 
code (IValenti &: Piskunov Ill996l ) and line list (IKupka et al.lll999l )) produce indistinguishable 
results. 

Secondly, we have applied spectral disentangling. This technique solves simultaneously 
for the contributions of the components to the observed composite spectra, and for the 
Doppler shifts in the component spectra. We have applied the method introduced by 
Simon fc SturmI (119941 ) using a version of the original code revised for the Linux operat- 
ing system. It assumes a constant light level outside eclipses, which is fulfilled for V12. 
Dedicated IDlfl programs were applied to remove cosmics ray events and other defects, and 
for normalization of the spectra before these analysis. Spectral disentangling is applicable 
because sufficient (16) WIYN spectra (S/N ^ 15-25) are available giving good phase/velocity 
coverage outside echpses. We have assumed the ephemeris given in Equation [2] (Section 13.31) 
and fitted for the velocity semi-amplitudes of the spectroscopic orbit for each component. 
We obtain Kp = 73.0 kms~^ and Kg = 74.7 kms~^, i.e. excellent agreement with the results 
given in Table [2l We have not quoted uncertainties for the semi- amplitudes since it is no t 
clear how to calculate robust values in the disentangling analysis (IHynes &: Maxtedl 119981 ). 
The disentangled spectra have S/N ratios of ~60 and cover 200A. 



''http: / /www. ittvis.com/ProductServices/IDL.aspx 
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Fig. 1. — Systematic errors in the raw TODCOR velocities of V 12 determined from simu- 
lations with synthetic binary spectra (filled circles: primary; open circles: secondary). The 
differences are plotted both as function of orbital phase (upper panel) and radial velocity 
relative to the center-of-mass velocity (lower panel), and have been applied to the measured 
velocities as corrections. Phase 0.0 corresponds to central primary eclipse. 



-8- 



Table 1: Radial velocities for VI 2. 



HJD Phase RVp RVs {O - C)p (O - C), 



- 2 400000 




kms 


-1 


kms ^ 


kms ^ 


kms ^ 


50331.86241 


0.72352 


32, 


.6 


-112.3 


0.4 


1.4 


50559.93366 


0.78823 


31 


.4 


-112.5 


0.3 


0.1 


50614.66644 


0.20309 


-110, 


.2 


32.3 


-0.1 


0.8 


50614.91145 


0.24076 


-112 


.2 


35.1 


0.9 


0.5 


50615.77338 


0.37328 


-92, 


.3 


13.2 


-0.0 


-0.2 


50615.91243 


0.39465 


-84, 


.6 


5.3 


0.4 


-0.6 


50780.65039 


0.72221 


31, 


.5 


-114.1 


-0.5 


-0.5 


50797.99374 


0.38866 


-86, 


.7 


8.4 


0.4 


0.3 


50815.71183 


0.11272 


-88, 


.7 


8.3 


-1.0 


-0.3 


50920.96276 


0.29447 


-111 


.0 


30.7 


-0.7 


-1.2 


50976.69696 


0.86330 


14, 


.8 


-97.0 


-0.6 


-0.4 


50996.77643 


0.95041 


-17, 


.1 


-63.5 


0.5 


-0.6 


51127.70306 


0.07966 


-75, 


.6 


-4.7 


-0.5 


-0.6 


51128.67943 


0.22977 


-112, 


.7 


35.1 


-0.1 


1.0 


51175.01367 


0.35341 


-97, 


.8 


19.3 


0.5 


-0.2 


52121.94531 


0.93896 


-11, 


.7 


-68.5 


0.9 


-0.5 
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Table 2: Spectroscopic orbital solution for V12. T is the time of central primary eclipse. 



Parameter 



Value 



Adjusted quantities: 
T (HJD-2 400 000) 
Kp (kms-i) 
K, (kms-i) 
7 (kms~-^) 

Adopted quantities: 

P (days) 

e 

Derived quantities: 

MsSin^i (Mo) 

Op sin i (10^ km) 
Os sini (10^ km) 
asini (R©) 

Other quantities 

Nobs 

Time span (days) 
CTp (kms~^) 
(Tg (kms~^) 



50906.0389 ± 0.0028 
73.2 ±0.2 
74.7 ±0.2 
-40.0 ±0.1 

6.5042969 


1.102 ±0.006 
1.080 ± 0.005 
0.980 ± 0.008 
6.547 ±0.019 

6.685 ±0.019 
19.011 ±0.040 

16 
1790 
0.65 
0.77 




Fig. 2. — Spectroscopic orbital solution for V 12 and radial velocities listed in Table [H (filled 
circle: primary; open circle: secondary). The horizontal dotted line (upper panel) represents 
the center-of-mass velocity of the system. Phase 0.0 corresponds to central primary eclipse. 
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2.3. Cluster Membership Probabilities for V 12 



The center of mass radial velocity of V 12 is —40.0 ± 0.1k m s~^ (Table [2]). From the 
woes radial-velocity survey of NGC 188, (iGeller et al.l l2008l ) report an average cluster 
velocity of —42.36 ± 0.04 kms~^. The center of mass velocity of V 12 is 2.5a from the fitted 
cluster mean, corresponding to a radial-velocity membership probability of 50%. A 97% 
as trometric membersh ip probability was determined for V 12 from the proper motion study 
of IPlatais et"aL] J2003h . 



3. Photometry of V 12 



The photometric data for V 12 consist of V and Ic CCD observations obtained at the 
2.56m Nordic Optical Telescope (NOT) July 2004 - April 2005, at the 1.2m Flemish Mercator 
Telescope April 2005, and at the 0.6m telescope at Bialkow observatory in Poland February 
2005 - April 2005. At the NOT, either StanCam or ALFOSC was used; we refer the reader 
to the telescope homepages for further information. Table [3] gives a summary of the data 
obtained. 

At each telescope the flat flelds used in the data reduction were obta ined during eve ning 
twilight. AU photometry wa s carried out with DAOPHOT/ALLSTAR dstetson lEoSTl ) and 
DAOGROW (IStetson Ill990l ) and transformed to a common coordinate system using MATCH 
and MASTER (P. Stetson, private comm.). For each frame we produced a point-spread 
function (PSF) using the brightest stars in the fleld and subsequently carried out aperture 
photometry in large apertures (neighbor stars subtracted) using the NEDA routine provided 
with DAOPHOT. Subsequently the aperture photometry was fed into DAOGROW to obtain 
the flnal large-aperture magnitudes. We found that this procedure gave slightly better 
photometric precision than using profile-fitting photometry. 



Table 3: Summary of photometric data obtained for V12. 



Telescope 


Nv 


Nj 




^expi^-^) 


NOT/StanCam 


131 


118 


120 


90 


NOT/ALFOSC 


197 


189 


60 


45 


Mercator 





40 




180 


Bialkow 


51 


36 


300 


200 
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3.1. Transformation to V and Ic 



For the analysis of tlie liglit curve it is important that all data in a given band are 
on the same photom etric system. For this purpose we used the photometry provided by 
Stetson et al. I (120041 ) which gives calibrated UBVRI photometry for a large fraction of the 
stars in NGC 188. We no te that R and / are on the Cousins system. The stars in common 
with lStetson et al. I (12004 ) were then used as local standards and the following transformation 
equations was adopted (for both V and J): 



mobs = Mstd + a X {B - /)std + P 



Capital letters denote calibrated magnitudes and lower case refers to instrumental mag- 
nitudes. For each combination of telescope, instrument, and filter, the values of a and /3 were 
found by the following procedure: 1) We fitted eq. [T]to the data from each frame (excluding 
V12 which is variable); 2) We calculated the average value of a, as it is expected to be 
constant in time; 3) Using the average value for a the values of jS were re-determined as the 
final frame zeropoint (/5). 

Using these values we calculated V and Iq magnitudes for V12. In the calculation of 
the zeropoints stars with [B — I) in the interval between 1.2 and 2.6 were used. Most of the 
stars had a color within ±0.3mag of V12. Thus the zeropoints are essentially based on stars 
of very similar color to V 12, ie. stars similar to the turnoff color of NGC 188. 



Before combining the photometry and adjusting zeropoints (see Section ing , we formed 
mean values of the out-of-eclipse observations of V 12 from ALFOSC a nd StanCam in order 
to obtain a higher precision than obtained by IStetson et al.l (120041 ). This is important 



for the determination of effective temperatures and distances. RMS errors of about 0.00 5 
magnitudes were obtained, compared to 0.013 (V) and 0.018 (/) for IStetson et al. I (120041 ). 
The mean V magnitudes from the two instruments agreed, whereas a difference of about 
0.02 magnitudes was seen in /. We adopt the following out-of-eclipse results for V12: V = 
14.745 ± 0.005, I = 13.951 ± 0.010, V -I = 0.794 ± 0.011. 



3.2. V12 light curves 

When combining calibrated photometry from several detector /filter combinations the 
photometric scale and zeropoints are always a concern. We therefore compared the phased 
light curves to check their quality. We therefore calculated the phased light curve to check its 
quality. The V photometry showed good consistency between the various sources with typical 
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offsets from the average of ±0™01. However tlie situation for the / filter was somewhat poorer 
with ALFOSC and StanCam data showing a consistent offset of (K°025mag. We have not been 



able to track down a definitive cause for these inconsistencies but note that Stetson et al. 



(120031 ) carried out an extensive analysis of a large dataset for the old open cluster NGC 
6791 and found that combining calibrated photometry from different telescope/filter/CCD 
combinations and different nights can exhibit scatter of the order ±0?^02. We note that the 
/ flat field of ALFOSC contains relatively large gradients near the center of the field where 
V 12 is located in our images. This instrument is a combined focal-reducer and spectrograph 
and the gradient seen is probably due to sky concentration. We speculate that this is related 
to the offsets between the ALFOSC and StanCam data. Offsets of similar size were found 
for the Mercator and Bialkow data in I. 

Since the internal precision of our data is significantly better than 0™01 we decided to 
adjust the zeropoints on a nightly basis for both the V and I data in order to produce a 
phased light curve with the smallest amount of s catter. To do this we first decided to adopt 



the values of V and I from IStetson et al.l (l2004j ) as our out-of-eclipse magnitudes for V 12 



since this is well calibrated and multi-epoch photometry (at random phases for V12). For 
each night wi th photometric data out (or partly out) of eclipse we adjusted the zeropoint to 



the values of IStetson et al. I (120041 ). Using the phased data, small offsets were then applied 



to the remaining light curve segments. 

It is important to emphasize that the depths of the eclipses of the final light curves are 
the same as observed with the different instruments. The data from Bialkow, NOT/StanCam 
and NOT/ALFOSC cover parts of the eclipse which include the central eclipse and out-of- 
eclipse data, and the observed eclipse depths in these agree to 0™01 to 0™02. We also note 
that for the analysis of the light curve the absolute level of the out-of-eclipse part does not 
have any influence. This only enters into the determination of the absolute magnitude of 
the stars (or equivalently the cluster distance modulus). The out-of-eclipse data from NOT 
show a scatter around a constant level of 3.4 mmag for both the V and I filter. 

The light curves are shown in Figs. [3] and HI A total of 379 and 383 measurements were 
obtained in V and /, respectively 0. Both eclipses have been covered on several nights, and 
there are enough data points outside eclipses to define the maximum light level because V 12 
is so well detached that essentially no out-of-eclipse variations are present. The eclipses are 
narrow and of almost identical depths (about 0.6 mag in V). They are of the same duration, 
and secondary eclipse occurs at phase 0.5, consistent with a circular orbit. 

As a final note we point out that the V 12 field is uncrowded and the stars are bright. 



Light curve tables are available electronically 
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This means that it should be a fairly easy task to obtain even better data than presented 
here. To avoid the transformation difficulties encountered here this should be done using a 
single telescope/detector/filter combination for each band. 



3.3. Times of minimum and ephemeris 

From the V and I light-curve observations, six times of primary minimum and four 
times of secondary minimum have been derived. They are given in Table H] together with 
the time of primary minimum obtained from th e spectroscop i c orb it (see Section [2]) and 



the appro ximate time of seconda r y min imum by IZhang et al. I (120041 ). When possible, the 



method of iKwee fc van WoerdenI (119561 ) was applied, otherwise a second order polynomial 
was fitted to the observations. 

An unweighted linear least squares fit to the times of primary minimum yields the 
following ephemeris: 



Min I = HJD 2453299.62009 

±25 



+ 6'?5042969 x E 
±18 



(2) 



which we adopt for the light-curve and radial velocity analysis. Within errors, the times 
of secondary minimum give the same period, and applying weights leads to nearly identical 
results. We find no evidence for eccentricity of the orbit of V12. 



3.4. Photometric elements 



The choice of binary model for light-curve analysis depends strongly on the relative radii 
of the components and the irradiation between them. In our case, the system is well detached 
with relative radii below 0.1, the components are little deformed (^oblatenes s of 0.0006) and 
refiection effects are small. Theref ore, the simple Nelson-Davis-Etzel model (INelson fc Davis 
19721 : lEtzellll98ll : iMartynov Ill973l ). which represents the deformed stars as biaxial ellipsoids 
and applies a simple bolometric refiection model, is expected to be adequate. 



We have adopted the corresponding EBOP code (jPopper fc Etzellll98ll ). supplemented 
by an extended version JKTEBOP including Monte Carlo simulations fo r the analysis and 



the ass ignment of realistic errors for the photometric elements; see e.g. ISouthworth et al. 
(l2004al lbl). 



We have analyzed the V and / light curves independently, and for each band identical 
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Fig. 3. — Observed V light curve of V12 (upper panel, crosses) and the theoretical EBOP 
light curve (line) corresponding to the photomet r ic sol ution presented in Table [5] assuming 
linear limb darkening coefficients by I Van Hamme Jl993h . (0-C) residuals of the observations 
from the theoretical light curve are shown in the two lower panels. 
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Fig. 4. — Observed I light curve of V 12 (upper panel, crosses) and the theoretical EBOP 
light curve (line) corresponding to the photomet r ic sol ution presented in Table [5] assuming 
linear limb darkening coefficients by I Van Hamme Jl993h . (0-C) residuals of the observations 
from the theoretical light curve are shown in the two lower panels. 
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Fig. 5. — E BOP solutions for V 12 for a range of fixed k values. Linear limb-darkening 
coefficients by IVan Hammd (119931 ) were adopted. The upper left panel shows RMS errors of 
the fit to the observations. Symbols are: cross V; square I. 
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Fig. 6. — Best fitting parameter values for tlie 10 000 syntlietic V 12 / ligfit cu r ves cr eated 
for the Monte Carlo analysis. Linear limb darkening coefficients by I Van Hammd (119931 ) were 
adopted. 
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weights have been assigned to all observations. A photometric scale factor (the magnitude 
at quadrature) was always included as an adjustable parameter, and the phase of primary 
eclipse was allowed to shift from 0.0. The mass ratio between the components was kept at the 
spectroscopic value. Gravity darkening coefficients corresponding to convective atmospheres 
were applied in accordance with the effective temperatures derived below, and the simple 
bolometric reflection model build into EBOP wa s used. Linear lirab darkening coefficients 
were either assigned from theoretical calculations (jVan Hammdl 19931 : IClaret Il2000l ) according 
to the effective temperatures, surface gravities, and abundances, or included as adjustable 
parameters. 

In tables and text on photometric solutions we use the following symbols: i orbital incli- 
nation; r relative radius; k = rg/rp] u linear limb darkening coefficient; y gravity darkening 
coefficient; J central surface brightness; L luminosity. 

EBOP solutions for V 12 are shown in Table [51 The following parameters were adjusted: 



I, r. 



k, and Js/Jp, and as seen their formal errors are quite small. For each band, the 
solutions for the different assumptions on limb darkening agree well, and also differences 
between the V a nd / solutions are sm all. The adjusted l inear limb dar kening coefficients are 
close to those by I Van Hammd (Il993l ). whereas those by I Claret I (120001 ) are about 0.1 higher 
in both bands. 

In order to assign realistic errors to the photometric elements, we have supplemented 
the analysis with EBOP solutions for a range of k = Vg/rp values near 1.00, which are 
compared in Fig. [5l and 10 000 JKTEBOP Monte Carlo simulations in each band; see Fig. [6] 
for details. The adopted photometric elements listed in Table |6] are the we ighted mean 
values of the EBOP solutions adopting the linear limb darkening coefficients by I Van Hamme 
(119931 ). The final radiative parameters were derived from a solution that uses the weigthed 
average of the geometric properties. Errors are based on the Monte Carlo simulations and 
comparison between the V and I solutions. We find that at phase 0.0, about 80% of the V 
light from the primary component is eclipsed, and at phase 0.5 about 85% of the V light of 
the secondary is blocked out. 



4. Physical parameters for V 12 

Absolute dimensions for the components of V 12 were calculated from the spectroscopic 
and photometric elements given in Tables [2l and [61 As seen in Table [71 masses and radii for 
the only slightly different components have been determined to an accuracy of about 0.6% 
and 1.4%, respectively. Individual V and / magnitudes are included, as calculated from the 
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combined magnitudes outside eclipses (Section 13. ip and the luminosity ratios (Table [6]). For 
interstellar reddening and absorption corrections, we have used E{B — V) = 0^087 ±0^010, 
E{V -I) = 1.3 1 X E(B - V) = ^l] 4:, and Ay = 3.1 x E{B -V) = 0""27. For Ic (8140A), 
the relations by lCardelli et al. I (119891 ) give the same EjV — I)/E { B — V ) ratio, and for the 
color of V12, the Eiy - I)/E{B - V) relation by IPean et all Jl978h leads to the same 
{V — I) reddening. 



4.1. The systemic and component effective temperatures 



Ef fective te r apera t ures Tf>ff were derived froni the recent calibrations by [A 



( 1996 ). ICasteUi lfll999h.lVandenBerg fc ClemI (120031 ) . iRamfrez fc Melendez I (l2005f ). lMasana et al 



onso et al. 



( l2006l ). and ICasagrande et al. I (j2006l ). Because the surface fluxes of the components are 
nearly i dentical, the tw o stars have almost the same T^ ff; the empiric al visual flux calibra- 



tion by iPopper I (jl980l ) with bolometric corrections by iFlower I (119961 ) yield a difference of 
only 25 K. We have therefore ap plied all available color indices for the c ombined light from 
Section l3Jllstetson et all tooi ). and from 2MAsil JHKs photometry JCutri et al.ll2003h . 
and the individual {V — I)o indices listed in Table [3 Results from the combined indices are 
presented in Table [H For each color index, the temperatures derived from the different cali- 
brations agree quite well. Differences in the derived temperatures between indices is, in part, 
due to uncertainties of the indices, which cause changes of 50-100 K, except for ~250 K in 
{R — I). The temperatures based on (V — R) and (R — I) would increase by about 120 K and 
decrease by about 20 K, respectively, if the possible zero point error in R of ™0.02 (noted by 
Stetson et al. I (120041 )) is applied. A somewhat better agreement between the temperatures 
is then obtained from the different indices. 

For [Fe/H] = 0.( 30, the individual (V — I)n in dices hsted in Table [7] give temperatures of 
5818 K and 5750 K jRamirez fc Melendez Il2005h 5928 K and 5858 K lvandenBerg fc Cle^ 
20031 ) . and 5829 K and 5765 K (ICasagrande et al. 1120061 ) for the primary and secondary com- 
ponents, respectively. Formally we obtain temperature differences between the components 
higher than 25 K, as derived from the V surface flux ratio, but the uncertainties of the 
indices, the reddening, the metal abundance, and the calibrations themselves add up (in 
quadrature) to individual temperature uncertainties of about ±100 K. Using [Fe/H] = -0.1 
typically decr eases the temperatures by l ess than 100 degrees. Giving higher weight to the 
calibration by I VandenBerg fc ClemI (120031 ). we assign 5900 K and 5875 K to the primary and 
secondary components, respectively, both with uncertainties of ±100 K. 
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4.2. Analysis of disentangled spectra 



Analysis of the disentangled spectra was performed to provide independent values for the 
temperatures, metallicities, and rotational velocities of the components. We cross-correlated 
the disentangled spectra of each binary component against a grid of synthetic spectra0 span- 
ning appropriate ranges in effective temperature, surface gravity, metallicity, and rotation 
velocity. The best matches (highest cross-correlation peak heights) for fixed surface gravities 
(logg) of 4.17 and 4.19 were obtained for temperatures of 5900 K and 6000 K, metallicities 
of -0.14 and -0.13, and rotation velocities of 15kms~^ and 17kms~^, respectively, for the 
V12 A and B components. The same analysis was performed for 5 twilight sky spectra 
from WIYN/Hydra, resulting in a mean temperature 60 K above solar, a mean metallicity 
0.06 below solar, and a mean rotational velocity 3.5km s~^ above solar. We estimate that 
the uncertainties on temperatures and metallicities derived from WIYN/Hydra spectra from 
this technique are at least 150 K and 0.1 dex. Thus considering the systematic offsets for the 
Sun , the spectral analysis of V 12 gives results that are entirely consistent with the results 
from the photometric analysis. Furthermore, the rotational velocities of the binary compo- 
nents are consistent with the synchronous rotation (11km s^^) expected in an old circularized 
binary system. 



4.3. The distance and age of V 12 



The accurate masses and radii determined for the components of V 12 are independent 
of distance, reddening, and color-temperature transformations, and allow direct comparison 
with stellar models. From the absolute dimensions of the components, we determine identical 
distance moduh of Vo - My = ll'^2A ± 0^09 ((m - M)v = 11"'51 ± 0°'09), corres ponding 



to 177 ± 75 pc. We adopt uncertainties of ±0.02 for the bolome tric corrections by [Flower 



( 1l996l ). From main-sequence fi tting in the 
adopting E{B - V) 



Stetson et a 



t004 ] {B 



y)o. My CMD and 

0°'087, IVandenBerg fc StetsonI (120041 ) derive {m - M)v values of 



the 


Saraiedini f 


it al. 


Michaud et al. 


(2004) 



Based on 
= 0^^087, 



11™40 from models of solar metallicity. 



To determine the age of V 12, we have chosen to use the recent Y"^ models by lDemarque et al 



( I2OO4J ) and the latest VRSS Victoria-Regina models for scaled solar abundances for the 



^The library is based on ATLAS9 (Ihttp:/ /kurucz. harvard.edu/') and the companion program SYNTHE 
used to compute a synthetic spectrum from the model atmosphere and line list. The R ~ 500, 000 synthetic 
spectra were broadened to an instrumental profile equal to the WIYN spectra. 
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heavy elements (IVandenBerg et al.ll2006l ). Comparisons were done for [Fe/H] = 0.00 ((Y,Z) 
= (0.266,0.018)) and [Fe/H] = -0.10 ((Y,Z) = (0.259,0.014)), corresponding to the range 
of recent values for the cluster metallicity. The two sets of models treat convective core 
overshooting differently, but at the age of NGC 188 overshooting is not required to match 
the turnoff region w here V12 is located. We refer the reader to the original papers or 



Clausen et al. I (120081 ) for more details on the models and their input physics. 



Figure [7] shows the mass-radius diagrams with the locations of the V 12 components 
and the and VRSS isochrones for relevant ages and [Fe/H] of 0.0 and -0.1. Grids of 
isochrones with 0.1 Gyr age separations were computed for both models using the available 
interpolation routines. For [Fe/H] = 0.00, the 6.4 Gyr isochrones from both models, provide 
the best fit to the V12 components. For [Fe/H] = -0.1, the best-fit VRSS isochrone has an 
age of 6.0 Gyr, and the best fit Y"^ isochrone has an age of 5.9 Gyr. TableOlists the isochrone 
ages for the two different models and the two different chemical compositions considered. 
We adopt an age of 6.2 ±0.2 Gyr for V 12 and NGC 188. Thus, for solar metallicity, we agree 
with the age of 6.2 Gyr obtained by IVandenBerg fc StetsonI (12004r). We also c onfirm the age 
of 6.4 G yr derived from solar metallicity diffusive models by iMichaud et al. I (l2004l ) and the 



result by 



Salaris et al. 



mm which is 6.3 ±0.8 Gyr for [Fe/H] 



-0.03 ±0.06. In contrast, 
VandenBerg fc Stetson! (|2004| ) derive a much higher age (7.7 Gyr) than us for [Fe/H] = -0.10. 



Overall, the isochrone method applied to the NGC 188 CMD results in a much larger age 
uncertainty for V 12 than found by our direct method over the same abundance interval. 

In Figure [S] we present mass-luminosity comparisons with models. On the theoretical 
side, the luminosity is independent of the mixing length assumed for convection; this is not 
the case for the radius. On the observational side, the luminosity includes uncertainties 
for the temperature calibration. In general, the models that best fit the components in 
the mass-radius plane are too luminous, with only the solar metallicity 6.4 Gyr model 
marginally within the component errors. The best fit of VRSS and Y^ models to the masses 
and luminosities of the binary components gives ages of 6.1 Gyr and 6.3 Gyr for [Fe/H] = 
0.0, and 5.3 Gyr and 5.4 Gyr for [Fe/H] = -0.1. To reach agreement between the model ages 
determined from the mass-radius and mass-luminosity diagrams would require an increase in 
the effective temperature of the V 12 components. The necessary temperature changes are 
~150K and ~400K for the [Fe/H] = 0.0 and [Fe/H] = -0.1 VRSS isochrones, and ~760K 
for the [Fe/H] = -0.1 Y^ isochrone. We note that the high sensitivity of luminosity to 
the effective temperature greatly enhances the uncertainty on the luminosity, making ages 
derived from the mass-luminosity diagrams less reliable. 

Figure [9] presents the effective temperature vs. radius diagram comparing the values of 
the V 12 components to those of the VRSS and Y"^ isochrones and 1.1 mass tracks. In 
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Fig. 7. — a) Mass-radius diagram with V 12 components and best fit VRSS isochrones for 
[Fe/H] = 0.0 (solid line) and [Fe/H] = -0.1 (dotted line). The ages of the isochrones in Gyr 
are listed in the figure in the same left-to-right order, b) Mass-radius diagram with V 12 
components and best fit isochrones for [Fe/H] = 0.0 (solid line) and [Fe/H] = -0.1 (dotted 
line). Isochrone ages are listed as in a). 

Figure [9^ we show the 6.4 Gyr [Fe/H]=0.0 and the 6.0 Gyr [Fe/H]=-0.1 VRSS isochrones as 
solid and dotted curves, respectively. In Figure [9b the 6.4 Gyr [Fe/H]=0.0 and the 5.9 Gyr 
[Fe/H]=-0.1 isochrones are plotted as solid and dotted curves, respectively. In both 
figures mass tracks corresponding to 1.1 Mq are shown as grey solid and dashed curves for 
[Fe/H]=0.0 and [Fe/H]=-0.1, respectively. For solar metallicity, the VRSS isochrone and 
track only marginally fit the observed quantities of the V12 components, being slightly 
hotter, whereas for [Fe/H]=-0.10 both are clearly too hot, especially noting that abundance 
uncertainties of ±0.10 are included in the temperature errors. Similarly, the isochrone and 
track for [Fe/H]=-0.10 are also too hot for a proper match to the binary components, although 
slightly closer in temperature than the corresponding VRSS models. The solar metallicity 
Y"^ 6.4 Gyr isochrone and 1.1 mass-track, however, provide a very close match to the 
V12 components. We also note, that the point on the solar metallicity 6.4 Gyr isochrone 
corresponding to a stellar mass of 1.09 Mq, falls appropriately between the 1.103 Mq and 
1.081 components. 

From the mass vs. luminosity and effective temperature vs. radius comparisons, the 
6.4 Gyr solar-metallicity Y'^ model provide the best match to the observed and derived prop- 
erties of the V12 components. However, we discuss in Section O how a zeropoint error in the 
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Fig. 8. — Mass-luminosity diagram with V12 components. In a) the 6.4 Gyr [Fe/H] = 0.0 
and 6.0 Gyr [Fe/H] = -0.1 VRSS isochrones are shown as sohd and dotted curves, respectively. 
In b) the 6.4 Gyr [Fe/H] = 0.0 and 5.9 Gyr [Fe/H] = -0.1 Y'^ isochrones are shown as solid and 
dotted curves, respectively. The isochrone ages are those determined from the mass-radius 
diagram. 



I-band photometry from lStetson et al.l (120041 ). may cause the derived temperatures for the 
V12 components to be too cool by ~100-150K. With a correction for such a temperature 
offset, the VRSS and Y"^ isochrones and mass-tracks would fit V12 equally well. 



5. NGC 188 and V 12 



The purpose of this section is to present V12 in the NGC 188 CMD, and to investigate 
how well the VRSS and Y^ isochrones that best fit the V 12 components in the mass-radius 
diagram, fit the rest of the cluster's main- sequence, turn-off, and giant bra nch stars. W e 



show in Figure [TD] the NGC 188 V-I vs. V CMD using the photometry from IStetson et al. 



(120041 ) ■ All stars brighter than My o f 4.5 are single radia l- velocity and proper- motion mem- 



bers of NGC 188 JCeller et al.lbood : IPlatais et aPboO sh. Stars fainter than My of 4.5 



are 



proper- motion members of NGC 188 (IPlatais et al.l 12003 ). The position of V12 and of its 
components are marked using the V magnitudes and V-I colors listed in Table 7. Adopting 
the distance (Vq — My = 11™24) and ages of V 12 determined from the analysis of V12 in 
Section B73l and a reddening of 0.087, we overplot in the left panel of Figure [TOl the 6.4 Gyr 
solar metallicity and the 6.0 Gyr [Fe/H]=-0.1 VRSS isochrones. In the right panel of Fig- 
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Fig. 9. — The effective temperature vs. radius diagram comparing the values of the V 12 
components to those of VRSS and isochrones and mass-tracks, a) The 6.4 Gyr [Fe/H]=0.0 
and 6.0 Gyr [Fe/H]=-0.1 VRSS isochrones are shown as sohd and dashed curves, respec- 
tively. The solid and dashed grey lines represent the VRSS 1.1 evolutionary tracks for 
[Fe/H]=0.0 and [Fe/H]=-0.1, respectively, b) Black solid and dashed curves show the 6.4 Gyr 
[Fe/H]=0.0 and 5.9 Gyr [Fe/H]=-0.1 Y"^ isochrones, respectively. The Y^ 1.1 Mq mass-tracks 
are shown as in a). In both a) and b) the horizontal dashed lines mark the points on the 
solar-metallicity isochrones corresponding to a stellar mass of 1.09 M©. 



ure [To] we have overplotted the 6.4 Gyr solar metallicity and the 5.9 Gyr [Fe/H]=-0.1 Y^ 
isochrones. For the solar metallicity isochrones we have marked the turnoff masses nearest 
the V 12 primary and secondary components. 

Both VRSS isochrones fit the main- sequence, turnoff, and giant branch of NGC 188 
well. The Y^ isochrones follow the turnoff curvature well, but both are too blue on the 
lower main-sequence, and the metal-poor Y"^ isochrone also falls to the blue of the giant 
branch. It is clear that without the constraints on the distance and age set by V12, and 
on the cluster metallicity by spectroscopic studies, good fits to the cluster stars could be 
obtained from a range of di fferent cluster metalh cities, distances, and ages. For NGC 188, 



this is demonstrated well by IStetson et al. I (120041 ). who show that ages between 5.9 Gyr and 



8.1 Gyr result from the isochrone method when plausible ranges in metallicity and distance 
modulus are considered. This age-range corresponds to ~33% of their mean cluster age of 
6.8^g;9Gyr. With V12 we provide a constraint on the age of NGC 188 with only a <10% 
uncertainty due primarily to the uncertainty in the cluster metallicity. Judging from current 
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Fig. 10. — The color-magnitude diagram for NGC 188 based on photometry from 
Stetson et al. N2004h . All star s brighter than My of 4.5 are single rad ial- velocity and proper- 
motion members of NGC 188 JOeller et al. Iboosl : lpiat_ais et al.lbnoal) . Stars fainter than My 
of 4.5 are proper-motion members of NGC 188 (jPlatais et al. 1120031 ). The position of V 12 
and the individual components are marked using the colors and magnitudes listed in Table 
7. In the left panel we show the 6.4 Gyr [Fe/H] = 0.0 (solid) and the 6.0 Gyr [Fe/H]=-0.1 
(dashed) VRSS isochrones. Relevant mass-points are marked on the 6.4 Gyr ([Fe/H] = 0.0) 
isochrone. In the right panel, we display in the same manner the 6.4 Gyr [Fe/H] = 0.0 (solid) 
and the 5.9 Gyr [Fe/H]=-0.1 (dashed) isochrones. 
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spectroscopic measurements and from the best fit isochrones to the cluster, it is hkely that 
the cluster metallicity is between solar and = -0.1. With better determination of the cluster 
metallicity, the age of V 12 and NGC 188 can be further improved. 



We note that in the V-I colors from lStetson et al. I (120041 ). V 12 and thus both its compo- 
nents are offset to the red by a few hundreds of a magnitude relative to cluster turnoff. The 
offset is likely due to a zeropoint error in the I band, because the Stetson et al. B-V color of 
V12 and its components match the B-V color of the cluster t urnoff. Similar off s ets fo r V 12 
are found in the B-V and V-I colors in the photometry fro m Sarajedini et al. ( 19991 ). Be- 
cause we calibrate our photometry using the IStetson et al.l (120041 ) photometry for V 12 (see 
Section ISTTl) . we inherit this I band zeropoint shift. Therefore, in Figure [TOl V12 and both 
its components display a slight red offset with respect to the cluster turnoff. The V-I color 
offset correspond to ~100-150K and may explain, in part, the discrepancies in luminosity 
and effective temperature between V12A,B and the VRSS and Y2 isochrones in Figure [8] 
and Figure [H The color offset has no affect on the age determination for V12 in this paper 
because the age is based on the stellar masses and radii only. However, any change in the 
color of V12 will affect the distance determined for the binary and thus NGC 188. 



6. Summary and conclusions 

We have presented time-series photometric and spectroscopic observations of the de- 
tached, double-lined, eclipsing spectroscopic binary V12 in the old open cluster NGC 188. 
We have determined masses of 1.103 ± 0.007 Mq and 1.081 ± 0.007 Mq, and radii of 
1.424 ± 0.019 Rq and 1.373 ± 0.019 Rq for the binary primary and secondary components. 
The uncertainties in the component masses and radii are not affected by the reddening, 
metallicity, and the distance to the binary. From the absolute dimensions of the binary com- 
ponents, we derive identical distance moduli of Vq — My = 11™24±0?^09 for both components, 
corresponding to 1770 ± 75 pc. 

By comparing the masses and radii of the binary components to the Victoria-Regina 
(VRSS) and the Yale (F^) model isochrones with [Fe/H] = -0.1 and [Fe/H] = 0.0, we de- 
termine a mean age of 6.2 Gyr for V12 and thus NGC 188, with an uncertainty of 0.2 Gyr. 
This uncertainty is primarily due to the 0.1 range in [Fe/H] for NGC 188 found in t he lit- 
erature. Our age for V12 agrees with the age determined by VandenBerg fc StetsonI (2004) 



for a solar metall i city i sochrone, and with the the ages derived by iMichaud et al.l (120041 ) 



and ISalaris et al. I (120041 ). We find good agreement between the ages derived from both the 
VRSS and Y'^ models. The solar-metallicity 6.4 Gyr isochrone match well the masses, 
radii, and effective temperatures of both the primary and secondary components of V 12. 
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The corresponding VRSS isochrone have temperatures shghtly hotter (~ 30 K) at the com- 
ponent radii, but still fit the compone nts within the errors. A systematic offset in the V-I 
color of V12 from lStetson et al.l (120041 ) may explain the temperature discrepancies between 
the V12 components and the VRSS isochrones and mass-tracks. Any future correction of 
the V-I color for V 12 will change its distance. 

We use the distance and ages for V 12, together with best estimates for the metallicity 
and reddening of NGC 188, to investigate the locations of the corresponding VRSS and Y"^ 
isochrones relative to cluster members in the CMD. We find that those isochrones fit the 
cluster cluster giant branch, turnoff, and upper main-sequence reasonably well. However, 
equally good or better fits to the cluster stars could be obtained by adopting di f ferent model 
metallicities, distances, and ages. This was demonstrated well in lStetson et al.l (120041 ) where 
plausible changes in the model metallicity and distance result in a range of ages from 5.9 Gyr 
to 8.1 Gyr for NGC 188 - more than 3 times that resulting from our analysis of V12. 

The age of V 12 and NGC 188 can be further improved from better determination of the 
metallicity of NGC 188 from spectroscopic studies of main sequence, turn off, and subgiants 
in NGC 188. Additional radial-velocity measurements, and an effort to obtain improved light 
curves in more bandpasses, could also provide masses and radii with smaller errors, and thus 
put better constraints on the stellar models. Finally, additional precise eclipse photometry 
in e.g. the K-band can make possibl e an independent determination of the d istance to V12 



via surface-brightness relations (e.g. Idi Benedetto I Il998l : iKervella et al.ll2004l ) 
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Table 4: Times of minima for V12. 0-C values for primary (P) and phases for secondary (S) 
minima are calculated for the ephemeris given in Eq. O 



HJD-2400000 


r.m.s. 


Type 


0-C/Phase 


Band 


50906.03886 


0.00276 




0.00003 


Spectroscopic orbit 


53221.56770 


0.00052 


P 


-0.00083 


I 


53221.56854 


0.00053 


P 


0.00001 


V 


53229.61958 


0.00022 


P 


-0.00051 


I 


53229.62044 


0.00011 


P 


0.00035 


V 


53312.62875 


0.00019 


P 


0.00007 


I 


53312.62954 


0.00015 


P 


0.00086 


V 


51605.25 


0.01 


S 


0.49989 


Zhang et al (2001) 


53452.47085 


0.00028 


s 


0.49997 


I 


53452.47109 


0.00010 


s 


0.50000 


V 


53465.48005 


0.00047 


s 


0.50006 


I 


53465.48001 


0.00060 


s 


0.50005 


V 
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Table 5: Photometric solutions for V12 from the EBOP code. A photometric scale factor (the 
magnitude at quadrature) and the phase of primary eclipse were included as free parameters. 
Linear hmb darkening were kept at theoretical values or left free. COO = Claret (2000), 
VH93 = Van Hamme (1993). The errors quoted for the adjusted parameters are the formal 
errors determined from the iterative least squares solution procedure. 



Band 


V 


V 


V 


I 


/ 


/ 


Limb darkening 


COO 


VH93 


Free 


COO 


VH93 


Free 


I (°) 


88.55 


88.62 


88.65 


88.57 


88.62 


88.64 




±1 


±1 


±4 


±1 


±1 


±3 




0.0756 


0.0755 


0.0755 


0.0742 


0.0742 


0.0743 




±7 


±7 


±8 


±8 


±10 


±11 


rs 


0.0724 


0.0719 


0.0716 


0.0731 


0.0726 


0.0723 


k — Tg/Vp 


0.958 


0.952 


0.949 


0.985 


0.978 


0.973 


±17 


±18 


±20 


±23 


±26 


±29 




0.1480 


0.1474 


0.1471 


0.1473 


0.1469 


0.1466 


Up Ug 


0.67 


0.59 


0.55 


0.52 


0.43 


0.39 








±4 






±4 


Vv^Vs 


0.38 


0.38 


0.38 


0.30 


0.30 


0.30 


Jsj Jp 


0.9823 


0.9824 


0.9824 


0.9954 


0.9957 


0.9960 


±28 


±28 


±28 


±26 


±26 


±26 


Ls/ Lp 


0.901 


0.890 


0.884 


0.967 


0.953 


0.943 


a (mag.) 


0.0058 


0.0057 


0.0057 


0.0055 


0.0054 


0.0054 



Table 6: Adopted photometric elements for V12. The individual flux and luminosity ratios 
are based on the mean stellar and orbital parameters. 



i 


88?62 ± 0?05 




Tp 


0.0749 ± 0.0010 




Ts 


0.0722 ± 0.0010 






V 


I 


Jsj Jp 


0.981 


0.996 




±15 


±14 


Ls/ Lp 


0.913 


0.927 




±36 


±50 
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Table 7: Astrophysical data for V12. Wc have assumed T^/fQ = 5780 K, B.C.q = —0.08, 
and MboiQ = 4.74. We have adopted E{B-V) = 0.087±0.010, E{V-I) = 1.31xE{B-V), 
and Av = 3.1 x E{B - V). 



Primary Secondary 

Absolute dimensions: 



M/Mq 1.103 ±0.007 1.081 ±0.007 

R/Rq 1.424 ±0.019 1.373 ±0.019 

log^f (cgs) 4.174 ±0.012 4.196 ±0.013 

Photometric data: 

V 15.449 ±0.021 15.548 ±0.023 

Vo 15.180 ±0.037 15.278 ±0.039 

/ 14.663 ± 0.030 14.746 ± 0.032 

(V-I) 0.786 ±0.014 0.803 ±0.014 

{V-I)o 0.672 ±0.019 0.689 ±0.019 

Teff 5900 ± 100 5875 ± 100 

Mboi 3.88 ± 0.08 3.98 ± 0.08 

log L/Lq 0.34 ± 0.03 0.30 ± 0.03 

B.C. -0.06 -0.06 

My 3.94 ± 0.08 4.04 ± 0.08 

{m-M)v 11.507 ±0.085 11.504 ±0.086 

Vo-Mv 11.237 ±0.090 11.235 ± 0.091 

Distance (pc) 1768 ± 74 1766 ± 74 
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Table 8: Effective temperatures for th e com bined light of y i2 deri ved from recent call 
brations. AAM 
RM 



(l2003h. 



Alonso et al.l Jl996h 



C = ICastelh I t99%. VC = IVandenBerg fc Clem 



Ramirez fc Melendez I (120051 ). MJR = iMasana et al.l (120061) . and C LF 



Casagrande et al.l (120061 ) . Extinction ratios from Table 1 in iRamirez &: Melendez I (120051 ) 



were adopted. De-reddened colour indices are given below the observed values. 



{B-V) 


(V-R) 


(V-I) 


(R-I) 


(V-J) 


(V-H) 


{V - K,) 


Calibration 


0.687 


0.427 


0.795 


0.367 


1.291 


1.660 


1.658 




±16 


±18 


±11 


±20 


±29 


±33 


±31 




0.600 


0.375 


0.682 


0.306 


1.103 


1.442 


1.423 




5837 














AAM 


5980 


5705 


5899 


6146 








C 


5887 


5650 


5890 










VC 


5846 


5616 


5780 


5986 


5871 


5690 


5883 


RM 














5927 


MJR 


5958 


5659 


5791 


6011 


5874 


5804 


5969 


CLF 



Table 9: Ages for V12 as determined from fitting isochrones to the binary components in 
the mass-radius plane. Two different models and two different chemical compositions are 
considered. Age uncertainties due to mass and radius errors are about ±0.25 Gyr. 

Model [Fe/H] Y Z Age (Gyr) 

VRSS (Too 0.2768 0.0188 q7a 

-0.10 0.2684 0.0150 6.0 

0.00 0.2662 0.0181 6.4 

-0.10 0.2589 0.0145 5.9 



